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Abstract. We propose a novel one-dimensional simple model without disorder exhibiting slow dynamics 
and aging at the zero temperature limit. This slow dynamics is due to entropic barriers. We exactly solve 
the statics of the model. We derive an evolution equation for the slow modes of the dynamics which are 
responsible for the aging. This equation is equivalent to a random walker on the energetic landscape. 
This latter elementary model can be solved analytically up to some basic approximations and is shown to 
present aging by itself, as well as a slow logarithmic relaxation of the energy: (e)(t) ~ l/ln(t) at large t. 

PACS. 02.50.Ey Stochastic processes - 05.70.Ln Nonequilibrium thermodynamics, irreversible processes 



1 Introduction 

Below their glassy transition temperature, glassy systems 
relax very slowly and remain out of equilibrium on exper- 
imental time scales. Despite intensive theoretical studies, 
the exact nature of this transition is not clearly under- 
stood yet, even if it is widely believed that it is of dynam- 
ical nature [Q. Glassy systems seem to be trapped in some 
metastable states whose numbers and life times increase 
dramatically with decreasing energy. Therefore there is 
an increasing demand for paradigmatic toy-models con- 
taining the elementary physical mechanisms responsible 
for glassiness and its experimental manifestations, such 
as the aging effect 2 . Among several other mechanisms 
for slow relaxation and aging, the focus has recently been 
brought into the existence of entropic barriers 0. This 
term designates boundaries between regions of the phase 
space where the system has very rare possibility to find a 
path to go from one region to the other. The system can be 
trapped in such a region without the necessity of energetic 
barriers. By analogy with Arrhenius law, the height AS of 
an entropic barrier is defined by r = exp(AS), where r is 
the time the system needs to go through the entropic bar- 
rier. Simple models 0] - such as the backgammon model 
jSJ or urn models jH] - have been developed which ex- 
hibit such entropic barriers. These models are of mean- 
field type such as most of models for glassy dynamics. As 
compared to these simple models for glassiness and en- 
tropic barriers, the interaction in the present toy-model is 
one-dimensional and therefore it can be legitimately con- 
sidered as more physical. 

On the other hand, our model is based on the well- 
characterized permutation group S n , and many calcula- 



tions can be carried out either exactly or after some basic 
assumptions. 



2 The model 

We consider a one-dimensional system of size n where the 
configurations are the n\ permutations on {1, . . . , n} of the 
symmetric group S n . We represent a configuration a G S n 
by a word a = a(l)a(2) . . . a(n) with a(i) <E {1, . . . , n}. 
The energy E of a configuration is defined by its number 
of descents. One says there is a descent between i and i + 1 
if a(i) > a(i + l). For example, let n — 6, the permutation 



a ex = 15|236|4 



(1) 

and its energy is 



has 2 descents symbolized by an arrow 
therefore E(a ex ) = 2. 

We denote by D k the degeneracy of the energy level 
E = k. The D„ are known as Euler numbers[Jj. The iden- 
tity is the unique ground state with zero energy so that 

= 1. Furthermore we can define a symmetric permu- 
tation er c for all a € S n , by 



= n + 1 — cr(i), 



(2) 



for 1 < i < n, which has the energy E(a c ) = n — 1 — E(a). 
The distribution of configurations with energy k is sym- 
metric with respect to (n — l)/2, and = D™~ 1 ~ k . 
At high temperatures all configurations have the same 
probability, the mean energy per particle in this limit is 
(e) = (n — l)/2n, according to the above symmetry and 
(e) = 1/2 at the large size limit. 
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The system evolves via a Metropolis Monte-Carlo al- 
gorithm [H] and we study two different dynamics: in the 
non-local one, at each time step, any two sites i ^ j are 
chosen at random and we try to transpose the elements 
a(i), a(j) with a certain probability Pji, which depends on 
the temperature T and on the energy variation AE if the 
transposition were executed: Pji — min(l, exp(— AE/T)). 
In the local dynamics we choose only one site i at random, 
and we try to transpose a(i) and a(i + 1), with the same 
transition probability as above. Sometimes in the following 
we will call particles the o~(i). Indeed, this model can be 
seen as a system of n distinguishable labelled particles on 
a one-dimensional lattice which tend to sort themselves. 

We anticipate on the following to emphasize the ori- 
gin of the entropic barriers in the non-local dynamics. 
If quenched from infinite to vanishing temperature the 
system evacuates its energy and relaxes through decreas- 
ing energy levels towards its equilibrium state. We shall 
demonstrate in the following that it encounters entropic 
barriers between any two successive energy levels whereas 
it relaxes rapidly inside each energy level. For example let 
us explain why there is an entropic barrier between the 
first excited states and the unique fundamental one. We 
need to know more precisely the values of D^. Starting 
from the expression 0: 



(k + l) n = J2 Dh r> 
3=0 



71 + k — j 



(3) 



we derive the value of D\ = 2™ — (n + 1) ~ 2™ for n ^s> 1. 
We get in the same way: 



D k n ~ (k+l) n , for fc/n«l. 



(4) 



We can easily recover the result for D\ by constructing 
the permutations with one descent as follows: we place 
the elements of the permutation, beginning by 1 and con- 
tinuing with 2 etc., randomly in two different sets, like 
if we were distributing distinguishable particles into two 
boxes. We can build 2™ permutations that have a descent 
at the boundary of the two sets, except if we have built 
the identity, which is obtained in n + 1 ways, moving this 
boundary between any two successive elements of the iden- 
tity. Subsequently = 2" - (n + 1). This point of view 
can be generalized to any energy k. The system is seen 
as (k + 1) ordered subsets separated by descents. In this 
point of view, the descents can be considered as domain 
walls between ordered subsets. 

This point of view also highlights the similitude be- 
tween our model and the backgammon [HE] one. It consists 
in N distinguishable particles placed in N boxes, where 
the energy is equal to minus the number of empty boxes. 
Our subsets play the role of the boxes in the backgammon 
model. The difference with our model is that the parti- 
cle interaction in the backgammon model is not localized: 
it is intrinsically mean-field. Note however that there ex- 
ists a one-dimensional generalization of the backgammon 
model PHI where particles are allowed to jump to nearest 
neighbors only. However, this generalization complicates 
the analysis of the model ^0] ■ 



We have seen above that there are about 2 n first ex- 
cited states and only one fundamental. In order to find this 
state, the system with energy equal to one, will wander in 
the phase space until it finds a state connected to the 
fundamental one. If the wandering is uniform in the first 
level, the probability that the system finds such a state, 
will be (n(n — l)/2)/2", where n(n — l)/2 is the number 
of path starting from the fundamental state in the non- 
local dynamics. Being in one of these states, it will have 
a probability 2/n(n — 1) to choose the direction towards 
the minimal energy. So, there will be about 2™ steps of 
the dynamics before the system finds the identity state, 
which corresponds to an entropic barrier of height of order 
n. Note that the previous argument also holds in the local 
dynamics: one obtains the same entropic barrier. Never- 
theless, the type of the dynamics will be of importance, 
since one can suppose an uniform wandering only for the 
non-local one for which we can prove the connectivity of 
the energy levels even at T — 0. Which means that there 
is always a path between any two configurations of one 
level that is inside it. The system searching the paths to 
decrease its energy will not have to pass through any en- 
ergy barrier. It is not the case for the local dynamics. But 
the entropic barriers, which depend mainly on the ratio of 
the level sizes, will be also present in the local dynamics, 
in conjunction with energetic barriers. In terms of permu- 
tation, passing through an entropic barrier corresponds to 
a complex rearrangement of particle. The order in which 
pairs of particles must be swapped to optimize their posi- 
tion without increasing the energy is quite constrained. 



3 Statics 

One of the advantages of the model is that its partition 
functions, canonical or grand canonical, have been exten- 
sively studied in the mathematical literature. Indeed, we 
can write the canonical partition function as follows: 

n-i / , x n-i , 

z=5>* eX p (-1) = = i(i-tr +i 5> B * h : 



fe=0 



fc=0 



fc>0 



(5) 

with t = exp(— 1/T) and where the last equality comes 
from combinatorial analysis More precisely, it is ob- 
tained inductively using the following relation on the E)\: 



(n + l-fc)£>„ +(fe + !)£>. 



(G) 



This former equation is obtained by building configura- 
tions with n + 1 elements by adding the (n + l) th particle 
in all the possible positions on a configuration with n el- 
ements. One can further remark that Eq. comes from 
the second term of the right hand-side of Eq. © . To ob- 
tain the mean values of the thermodynamical observables 
we compute: 



In Z n + 1 , , , 1 , x - ^ , „ i- i 
= — - ln(l-i) + -ln> fc 71 ** -1 . 

11 11 T) ' J 



(7) 



fc>0 



V. DESOUTTER, N. DESTAINVILLE: Slow dynamics in the one-dimensional 'descent model' 



3 



The sum in the right hand side of the former equation 
reads 

^fc^-^^exp-nMfc + l)- A) ( 8 ) 

fe>0 fc>0 ' 

= ^exp(-n/(fe,T)) 

= (2vrnT 2 ) 1/2 exp(l/T~n), 

by taking the continuous limit, and evaluating the so- 
obtained integral by a saddle-point argument for n ^> 1. 
This continuous limit is only valid if the width of the 
Gaussian approximate of exp (— nf(k, T)), is large as com- 
pared to the spacing between two consecutive energy lev- 
els, Ak = 1, that is to say if T 3> l/2y/n. In the ther- 
modynamic limit it will be always true for any finite tem- 
perature. Using Eq. J7J) and thermostatics identities the 
thermodynamical observables can be exactly computed. 
The free energy per site obeys 

/ = -T(ln(l-t)+ln(nT)) (9) 

and the mean energy per particle is given by 

(e) = T — t/(l - t). (10) 

It follows that the specific heat is 

C v = l + </(T(l-t)) 2 , (11) 

and the entropy per particle is given by 

s = ((e) ~f)/T = hi(nT)+ln(l-t)-t/T(l-t) + l. (12) 

Note that this entropy is not extensive. Let us remark the 
existence of an interesting low temperature domain de- 
fined typically by T < 0.1, where the first term of Eq. J7|) 
is negligible so that (e) ~ T . It delimits a low energy 
domain, where ~ (fc + 1)™, which results in cntropic 
barriers. We will see in the following part that the slow 
dynamics takes place in this region of the configuration 
space. These results show that there is no thermodynam- 
ical transition at any finite temperature. 



4 Dynamics 

We have defined two types of dynamics above, for the 
sake of simplicity, we focus on the non-local one, and we 
shall discuss the local one in the conclusion. As we have 
seen above, we study the dynamics of the system after a 
quench from high temperature, during which the system 
tries to decrease its energy to reach a low temperature 
equilibrium state. In the following, we will focus only on 
low temperatures such that the system has to encounter 
cntropic barriers during its relaxation process. 

In order to analyse the dynamics, we first map our 
model on a random walker on the energy levels. It is on 
this model, which is much simpler, that we get analytical 



results. In a second part, we check numerically by Monte- 
Carlo simulations that the dynamics of the descent model 
is quite well described by this random walker. In order 
to justify this mapping let us make some hypotheses. We 
suppose that the equilibration time in any energy level is 
short as compared to the time spent to go through an en- 
tropic barrier between two energy levels. That is to say we 
suppose that there is no entropic barrier inside an energy 
level. This hypothesis is reasonable since the energy levels 
are connected, and will be corroborated in the following 
by numerical results. We also suppose that the paths be- 
tween two energy levels are uniformly distributed inside 
these levels. 

Following these hypotheses we compute the probabil- 
ity that the system goes from one level to another using 
the mean features of the levels. We take the slow dynam- 
ics of the model as equivalent to a dynamics between the 
energy levels. In other words, we map the dynamics of 
our model on a Markovian random walker on the energy 
levels, whose features are precised in the following. Since 
we know the at low energy, and the number of paths 
starting from any configuration, we know the total num- 
ber of paths starting from one level, D^(n— l)n/2. We are 
searching the number of paths allowing the system to de- 
crease its energy from one level A;. It is also the number of 
paths starting from one level I < k and arriving in k. Since 
the energy variation cannot exceed 2, / E {k — 2, k — 1}. Let 
us introduce a parameter p(k) which represents the pro- 
portion of paths starting from the level k—1 and allowing 
the system to increase its energy by one unit. It certainly 
depends on the energy level. We investigated numerically 
the configuration space of the model, and we concluded 
that for — = e < 0.1, p(k) ~ Ke, with K a constant that 
seems to depend very weakly on n. For n = 10000, we 
have K ~ 1.99 and K = 1.97 for n = 12. For the links 
between the fundamental state and the first excited ones 
(k = 1), one can trivially show that K exactly equals 2. 
In the following we will consider K = 2. 

So, p£)^ _1 (n — l)n/2 is the number of paths starting 
from one level k—1 and increasing the energy by one unit, 
and pD^~ 1 /Dn the proportion of paths starting from one 
level k allowing the system to decrease its energy by one 
unit. Being equilibrated in the level k, the system has to 
find one of these paths to decrease its energy. These events 
are very rare because of the great ratio -D„/D„ 1 at low 
energy. Indeed, the typical time Tk,i, in Monte-Carlo steps, 
to go through the entropic barrier (of height AS) between 
the energy levels k and k — 1, satisfies: 

r k .i = exp(AS) = ~ - — - — = - exp(n/fe). 

pDl p V k J P 

This time gives us the probability rate by unit step that 
the random walker goes from k to k— 1, u^k-i = 1/ T M- 
In the same way we can obtain the typical time — 
exp(2/(fc/?i— l/n))/q that the system goes from a level k 
to k — 2, q being the proportion of paths starting from k 
and arriving in fc + 2. Simulations show that q is of order 1. 
We obtain uJk^k-2 = 1/^,2, but in the following we will 
always neglect it since — ^k^k-i < L 
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The master equation for the random walker at T = is 
given by: 

P(k,t + l/n) = P(k,t)+ 

p{k + 1) exp(-n/(* + 1))P(* + 1, t) 

-p(k)exp(-n/k)P(k,t), (14) 

where i is in Monte-Carlo unit (one unit is equal to n 
Monte-Carlo steps). By using the energy per particle e 
and the continuous limit in time and energy, this equation 
reads: 

= ^(p(e)P(e, t) exp(-l/e)). (15) 

Thus one obtains the evolution of the mean energy per 
particle with time (e)(i), using p(e) = 2e: 



To check ea.l|19|l. we made some numerical simulations on 
the random walker and on the original model. The inter- 
est of the random walker is that its simplicity allows us to 
compute numerically its exact dynamics for large n and t 
(up to t = 10 12 ) using its transition matrix. We find very 
good agreement between analytical and numerical calcu- 
lations. Monte-Carlo simulations on the descent model are 
also in very good agreement with analytical calculations 
at large t, as we can see in fig. [21 These results corroborate 
a posteriori the hypothesis that there are no entropic bar- 
riers inside the energy levels, since otherwise the random 
walker would be faster than the descent model. 

Let us remark that the only analytically unknown pa- 
rameter is pie). In fact, this parameter is irrelevant for 
the asymptotic dynamics while p is of polynomial type: 
p = Ke a . Indeed, using this equation for p, eq. (fTTj) reads: 



d(e) _ 8P{e,t) 



, e de 

dt J at 

/■(n-l)/n 

= -2/ eP(e,t)exp(-l/e)de. (16) 

Jo 

We calculate this integral by developing eP(e, t) exp(— 1/e) 
around (e): 



d(e) 
dt 



-2(e) exp(-l/(e)) + 0(Ae 2 exp(-l/(e))/(e) 3 ), 



where one keeps only the zeroth order in e 
shall demonstrate that Ae 2 = ((e — (e)) 2 ) 



(17) 

e), since we 
\e) 2 /2n. In- 
deed we can derive time evolution of Ae 2 from eq. 114f> pay- 
ing attention to keep terms of order 1/n since they cannot 
be neglected in this case. As above, discrete sums are eval- 
uated using an integral and developing eP(e, t) exp(— 1/e) 
around (e): 



dAe 2 
dt 



= -± e -WAe 2 + M e -V(e> _ I e - 2 /< e >. 



(18) 



The first term describes a rapid decay of Ae 2 towards its 
quasi-equilibrium value (e) 2 /2n. Hence Ae 2 will be sta- 
bilized around (e) 2 /2n, which is checked numerically for 
the descent model and the random walker. By integrat- 
ing eq. (|17f) we obtain the time the walker needs to go 
from high energy to the level e = (e), t = eexp(l/e)/2 + 
C(e 2 exp(l/e)) and: 



(e>(*) 



1 



ln(2t) + ln(ln(2i)) 



ln(t) 



(19) 



This dynamics can also be found in a variety of theoretical 
models 0], such as the backgammon model and the oscilla- 
tor model ^21 > as weu as m compaction of granular media 
(models and experiments |13U14U15| : see conclusion). This 
result means that the system will never reach the funda- 
mental state at T = 0, and will always stay out of equilib- 
rium despite the absence of energy barriers. Its relaxation 
time in the thermodynamical limit is diverging. In the 
long time regime the dynamics will be slower and slower. 



dt 



= -if(e) Q exp(-l/(e))+0(Z\e 2 exp(-l/(e))(e)"- 4 ), 

(20) 



and hence the mean energy per particle becomes: 



(c>(*) 



1 



1 



In(ift) + (2 - a) ln(bi(irt)) ln(f) 



(21) 



For a = 2, the evolution in l/ln(i) is even exact at all 
times in the thermodynamical limit, not only asymptoti- 
cally, which is well verified in numerical simulations. 

The variance Z\e 2 does not depend on a, since eq. (|18H 
becomes: 



^ = -2(e) Q - 2 e- 1 /^zi e 2 + ^ e -i/<e>. 
dt n 



(22) 



We still have the quasi-equilibrium value Ae 2 = (e) 2 /2n. 
We will see below that the shape of p does not change 
the aging in the energy-energy correlation function either. 
It is not surprising: what dominates the dynamics is the 
exponential decay of the random walker transition rate 
with the inverse energy. 

In the following we will consider again p — 2e, since it 
mimics the best the descent model. 

To study the dynamics at T > 0, we have to consider 
the possibility that the walker increases its energy. From 
the definitions of the dynamics and of p and q, it follows 
that LUk^k+i = _pexp(-/3) and ujk~>k+2 = gexp(-2/3), 
where (3 is the inverse temperature. At low temperature, 
<C 1, so we need not consider the pos- 
sibility that the system jumps two levels in one step. In 
the same way as at T = 0, the evolution of the random 
walker position follows a master equation, which in the 
continuous limit, leads to: 

= -2(e) ex P (-i y)+ 2(e) exp M) ^ 

+ G(Zi e 2 exp(-l/(e))/( e ) 4 ). 

We remark that d(e)/dt = for (e) = 1//3 = T, which is 
the equilibrium energy for the descent model for T« 1. 
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At high T the random walker does not represent the de- 
scent model anymore since it evolves in a region where 
£)„ ^ (fc+1)". Thus we focus on the low T region. One can 
see from eq. (j2HJ), that while the system is far from its equi- 
librium energy the possibility that it jumps to an higher 
level by thermal activation will be negligible in front of its 
possibility to decrease its energy: for (e) ^> 1/(3 we have 
exp(— -rly) 3> exp(— (3). Far from its equilibrium state, the 

system is not influenced by the temperature and its dy- 
namics is the same as at T — 0. This defines a typical 
time = exp((3)/(2(3) needed by the system to reach 
(e) = 1/(3 with the dynamics at T = 0. Near equilibrium, 
temperature effects appear and we obtain the dynamics 
by developing exp(— l/(e)) around 1/(3 in eq. J52J). We 
find an exponential relaxation of energy with the same 
characteristic time Tp (see fig. |21 inset). 

We define two types of correlation functions in order to 
characterize the dynamics. One is based on the matching 
of two permutations, in others words on the proportion 
of sites at time t which bear the same particles as at a 
time t w : C a (t,t w ) = Ei i,(i+t.),ui(t,„)/' 1 . where t w is 
the waiting time after the quench. At high temperature, 
it can be proven that C a (t ^ t w ,t w ) — 1/n, at lower 
temperature we numerically observe that it is also the 
case. This correlation function shows rapid decay for any 
t w , with no aging effects. The second correlation function 
we consider is the energy-energy correlation: 

n , ^ _ (e(t + t w )e(t w )}-(e(t + t w )}(e(t w )) 
MMtu) ~ Ae^) ' (24) 

We studied it numerically for both models at different t w . 
For the random walker, at T = one can see in fig. [3] that 
the system ages. The correlation function seems to tend 
towards the scaling law: C e (t,t w ) ~ t w /(t + t w ) at very 
large t w . 

At T > we expect the dynamics to be the same as 
at T = while the system is out-of-equilibrium. As long 
as t < Tp we find numerically the same plots as at T = 
whereas we find time translational invariance for t + t w > 
Tp. The same results are found for the descent model but 
with statistical noise since we cannot use the transition 
matrix. These results show the strong importance of the 
observable we focus on to observe aging |17| . 

Now we propose a simple argument leading to the fol- 
lowing law at large t and t w : 



Ce (^) tw ) 



t a 



In 2 f„ 



t + t w hf(t + t w ) 



(25) 



As illustrated in figure this law accounts quite well for 
our numerical observations as far as the random walker is 
concerned. In particular, the logarithmic corrections en- 
sure the collapse of the curves for different values of t w 
on a straight line of slope very close to -I. For the descent 
model, the statistical noise does not allow to distinguish 
if the logarithmic corrections improve the agreement be- 
tween theory and numerical experiments. 

Our argument is as follows: since the dynamics on en- 
ergy levels is Markovian, the future of the random walker 



beyond t w (i.e. t > 0) only depends on the energy distri- 
bution at time t w . For the descent model it will be only 
true at large times, when the entropic barriers make the 
dynamics on the energy levels Markovian. 

Now we make the following hypothesis: the evolution of 
the mean energy (e)(t + t w ) only depends on the mean en- 
ergy (e)(t w ) and not on the precise shape and width of the 
energy distribution at t w . This hypothesis is supported by 
numerical simulations with different initial energy distri- 
butions and it is corroborated by the excellent agreement 
between theoretical and numerical correlation functions. 
Furthermore, it amounts to neglecting the terms of higher 
order in eq. (|17|l . as we made before, which is exact in the 
thermodynamics limit. 

A consequence of this hypothesis is that any shift Se 
at t w of the mean energy has the same effect on (e)(t + t w ) 
as a time delay St at t w , since it only consists of a shift of 
the initial condition. For small Se. St is such that: 



Se d(e) 
~8i 



dt 



1 



tqt) ill tq i 



(26) 



This point is illustrated in figure ^ Now we still denote 
by (e)(t + t w ) the mean energy without any shift at t w , 
and by (e\Se, t w )(t + t w ) the mean energy after an energy 
shift Se at t w (conditional mean). Subsequently 



(e\5e,t w )(t + t w ) ~ (e)(t + t w + St) 



(27) 



(e)(t + t w )+St ^-(t + t w ) 



(e)(t + t w ) + 5e 



In 2 L 



t + t w ln 2 (t + t w )' 
Now the energy-energy correlation function is 
(e(t + t w )e(t w )) = ^ e * P ( e * \tw)P{e,t + t w \e* ,t w ) 

e* e 

where P{e* , t w ) is the probability that the energy per par- 
ticle is e* at time t w and P{e, t + t w \e* , t w ) is a conditional 
probability, and 

(e{t + t w )e(t w )) (29) 
= e*P{e\ t w )(e\5e = e* - (e) (t w ), t w ) 

e* 

= J2e*P(e*,t w ) x 



(e)(t + t w ) + (e* - (e){t w )) - — - — 7 



In 2 t v 



t + t w hi (t + C) 



(e)(t + t w )(e)(t w ) + 



In 2 t. 



t + t w \u\t + t w ) 



The last sum is the variance Ae 2 of the energy distribu- 
tion at time t w . Hence we get eq. (|25|) . Note that the same 
kind of argument can be used to derive the same energy- 
energy correlation function for the Backgammon model 
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and is also in excellent agreement with numerical simu- 
lations that we perfeormed independently. Note that our 
correlation function is not the same as the one calculated 
in reference |lfi| . which explains that is does not have the 
same expression. With p = Ke a , this law is unchanged as 
we saw before for (e) (t) . 

5 Conclusion 

We have focused so far on the non-local dynamics. Let 
us now discuss briefly the local one. Its interest lies in 
the fact that it is truly one-dimensional since it respects 
the one-dimensional character of the model. We shall see 
that even if its analysis is complicated by the existence of 
energetic barriers (which freeze the dynamics at T = 0), 
the qualitative conclusions are the same as for the non- 
local case. Indeed, one shows numerically that below an 
energy per particle e m (n) 0.1, which depends slowly on 
n, nearly all states are local potential minima. Therefore 
at T = 0, the system is always stuck in these minima, 
and at T > it has to pass over energy barriers at each 
energy level in order to lower its energy. However, these 
energy barriers are always of height 1, and one needs a 
time of order exp(/3) to pass them. Therefore, the times 
needed to pass over energetic barriers and entropic barri- 
ers have the same order of magnitude and the system is 
not substentially slowed by energetic barriers, at least at 
T > 0. More precisely, we have measured typical times 
r'(n, T) by evaluating the constant diffusion of particles 
at equilibrium, at small but finite temperatures and we 
conclude that r'(n, T) is of dominant order n 2 exp(/3). As 
compared to the non-local characteristic times, this dy- 
namics is n 2 times slower than the non-local one. Note 
that this prefactor n 2 is also present at high temperatures 
where it can be proven rigorously that r'(n, oo) ~ n 2 |18) . 
At high temperature, the n 2 term is certainly due to the 
diffusion of particles over the whole system. At smaller 
temperature, one can think that it results from the diffu- 
sion of descents acting as walls between ordered domains, 
which allows the system to explore energy levels. As far 
as aging is concerned, our numericals results remain com- 
patible with the law: C e (t, t w ) = t w j (t + t w ). 

A natural continuation of the present work will be to 
investigate into deaper detail the relationship between our 
model and compaction of granular media. Indeed, equa- 
tion H17[l also governs the evolution of the density in sim- 
ple models of compaction (see for example). On the 
other hand, in these mean-field models based on a free 
volume argument, as well as in the one- dimensional de- 
scent model, the slow dynamics is due to the necessity of 
complex and long rearrangements of particles to optimize 
the organization of the system. 

To finish with, we mention that as soon as a model 
presents an effective Markovian dynamics between energy 
levels due to entropic barriers, the present analysis can 
be applied to this model. For example, the Backgammon 
model can be tackled in such a way and one obtains the 
correct laws for the mean energy (e)(t) and the correlation 
function G e (t,t w ). 



We thank Remy Mosseri and Alexandre Lefevre for fruitful 
discussions. 
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Fig. 1. Schematic picture of the aging analysis. The full line represents the mean energy per particle whereas the dashed line 
represents the mean energy per particle when the system is at energy e* at time t w . St is the constant time shift between these 
two curves. 
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Fig. 2. (e)(t) for n — 500 and T = 0. Circles represent Monte-Carlo simulations of the descent model compared with the 
evolution given by l/(ln(Kt) + ln(ln(Kt))) (full line) which represents the energy of the random walker with p — Ke, K — 2. 
Circles in the inset show (e) for the descent model at a fixed time t — 10 4 = T(p=i2) for various temperatures and the full line 
the exact static energy. We check that for [3 such that Tp < 10 4 the system has reached its equilibrium energy, and that for all 
(3 such that > 10 4 the system is out-of-equilibrium and its energy is independent of (3. 
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Fig. 3. Energy-energy correlation of the random walker at T = 0, for n — 2000, as a function of (t + t w )/t w , compared with 
t w /(t + t w ) (full line). The inset shows the relaxation time as a function of /3, of the random walker, obtained from C e (t, t w ) at 
equilibrium, and the full line a fit by A exp(/3)/2/3 = Arp, with A — 1.25, corroborating analytical calculations. 
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Fig. 4. Energy-energy correlation corrected by In 2 (t + t w ) / In 2 (t w ) of the random walker at T = 0, for n = 2000, as a function 
of (t + t w )/t w , compared with t w /(t + t w ) (full line). The logarithmic corrections ensure the collapse of the data for different 
values of t w . 



